Two-spin entanglement induced by electron scattering in nanostructures 
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We present a model where two magnetic impurities in a discrete tight-binding ring become en- 
tangled because of scattering processes associated to the injection of a conduction electron. We 
introduce a weak-coupling approximation that allows us to solve the problem in a analytical way 
and compare the theory with the exact numerical results. We obtain the generation of entanglement 
both in a deterministic way and in a probabilistic one. The first case is intrinsically related to the 
structure of the two- impurity reduced density matrix, while the second one occurs when a projection 
on the electron state is performed. 
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I. INTRODUCTION 

The generation of entanglement in mesoscopic struc- 
tures is considered as a fundamental resource for the im- 
plementation of solid-state quantum information process- 
ing devices 1 . The first proposals for spin-based quantum 
computation concern direct interaction between qubits 2 . 
Further, different schemes for mesoscopic structures have 
been suggested that create separated streams of entan- 
gled particlea 3 i 4 i 5 i 6 i 7 i 8 . 

Recently, Costa et al*L have examined the possibility 
of entangling two spatially separated stationary spins by 
means of electron scattering. In this case the genera- 
tion of entanglement would require lower external con- 
trol. After considering a toy model concerning a bal- 
listic electron interacting in succession with two distant 
spins, Costa et al. face the more realistic problem of 
how two magnetic impurities embedded in a solid become 
entangled because of the injection of a conduction elec- 
tron, which is scattered by the impurities according to 
a s-d Hamiltonian. Since in a normal metal described 
by a tight-binding model the energy spectrum consti- 
tutes a continuous band, the appearance of an imaginary 
part in the eigenenergies of the system is expected that 
would limit the coherence time of the entangled state. A 
way to overcome this instability can be represented by 
the introduction of an artificial discrete system, i.e. a 
ring of iV-coupled quantum dots (QDs). During the last 
decade, electronic transport properties through quantum 
dots have been widely considered both experimentally 
and theoreticalljiiii. Because of the progress of nanotech- 
nology, it is possible to fabricate various structures of 
coupled QDs smaller than the electron coherence length. 

If the interaction term between the electron spin and 
the impurities is much less than the energy separation be- 
tween consecutive eigenvalues, in a finite-size system such 
as a nanostructure, resonance conditions are reached, and 
a reduction to a few-body systemAi can be observed. In 
that case an oscillatory regime is expected to come out 
also in the degree of entanglement. Dissipation effects 
could appear only through the interaction with some ex- 
ternal bath (for instance, coupling with phonons). 



The aim of this paper is to show how entanglement 
through electron scattering can be generated efficiently 
in such nanostructures. Then, we consider a finite tight- 
binding model, and add two magnetic impurities. By 
studying the dynamical evolution of the state of a con- 
duction electron injected in the chain, we establish the 
amount of entanglement between the two spins as a func- 
tion of time. We find that entanglement can be generated 
in a deterministic way as well as in a probabilistic one. 

The plan of the work is the following. After the In- 
troduction, in Sec. [H]we define the general Hamiltonian 
model, getting the equations of motion for the states in- 
volved in the evolution. In Sec. IIIII we establish the 
approximation of considering the interaction of the in- 
coming mode only with the other resonant modes, and 
solve the equation of motion. In Sec. IIVI the genera- 
tion of entanglement is studied by means of two different 
approaches: the measure of concurrence associated to 
the reduced density matrix obtained by tracing out the 
electron spin, and the "localizable entanglement" derived 
from a projective measurement on the electron. In Sec. 
El we conclude the paper. 



II. MODEL 

We consider a discrete ring structure, described by a 
standard tight-binding Hamiltonian, where two particu- 
lar sites (for instance, we label the first site with and 
the second one with L) are substituted by magnetic im- 
purities. The s-d Hamiltonian describing the system is 

H = H a + J - (§ ■ a Q + S L • , (1) 

with H = J2k,a e k a l,a a k,<7, where a\ a (a k ,a) creates 
(annihilates) one electron with spin a =|, j on the mode 
k, tk = —2wcosk are the eigenvalues of H in the ab- 
sence of spin interaction (w is the hopping amplitude 
between adjacent sites, k = (2n/N)n, N is the total 
number of sites, and n is an integer running from ~N/2 
to [(N/2) — 1]), J is the coupling constant between the 
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FIG. 1: Plot of the probabilities of the states involved in the 
evolution as a function of time. Solid lines correspond to theo- 
retical predictions, while dashed lines show numerical (exact) 
evolution. Red lines concern the state |TT) \[k), black lines are 
related to |||) ||_fc), while green lines regard |||) \\k)- The 
probabilities associated to |||) |T-0, IT-l) ltfc)> and IT-l) IT-fc) 
are not plotted, being practically indistinguishable from that 
°f I IT) Itfc)- The system parameters are the following: w — 1 
is chosen as unit of energy; the number of sites is N = 16; the 
scattering amplitude is J = 0.2, while the distance between 
the spins is L — 4. 
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where the notation |) refers to a configuration, while |) 
refers to the evolution of that configuration, and fi — 
1. The other states involved in the evolution satisfy the 
following equations: 
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impurity spins So and Sl and the electron spins uq and 
<?l, whose operators are defined as 



III. RESONANT COUPLING 
APPROXIMATION 
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Let us consider the impurities initially with spin up, and 
the introduction of one excess electron in the state ||fc). 
Then, differently from the treatment given in Ref. H we 
look at entanglement generation from a dynamical point 
of view, i.e., we analyze the evolution in the time do- 
main of the state |TT) life) to show explicitly how coher- 
ent effects persist. By introducing the complex Laplace 



The problem is significantly simplified by intro- 
ducing the following weak-coupling approximation. 
In the solution of the system derived from Eqs. 
©-0, it would appear denominators with the 

r i ~ 1 

structure ]u - e k - (,J 2 /N 2 ) J2 q f (*-«)/ (w - e q )\ , 

where / {k — q) is some weight function derived from 
(±1 ± e l ( 9 ~ fc ) L ) . If the scattering amplitude, of the order 
of J/N, is much less than the energy differences appear- 
ing in the spectrum of Hq, which is about w/N near the 
middle of the band, we keep just the resonant terms, cor- 
responding to q — ztk 11 . For instance, 
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and in the right-hand side of Eqs. ©~CJ we maintain 
only states with momentum k or —k, that is, the states 
that were degenerate with the initial one in absence of 
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interaction. Under these assumptions the equations of motion reduce to 
I 



imiu> 
imiT*) 



i 



W - £fe + jf 
1 



ITT) IU>- 2^(1 + e 



J 



^ - £fe + ^ ^ 



[lit) (IT* 



- 2lfcL )ITT)ll- fc ) 

)j + \n){\m 



-2ikL 



IT- fc )J] , 



i 



W - £fe 
1 

W - £fe 



liT) IT fe ) + ^ (-1 + e- 24fei ) liT) |T-fc) u + ^ ITT) (life), + li-fe) J 



ITDITa 



J 

27V 



(i 



J 



)IU)IT- fc ). + ^ITT)(IUL + e - 24fcL ll- fc )J 



(9) 
(10) 

(11) 



Together with these equations we must consider also instance, if k = n/2 and L is even, e ±2lkL = 1. In this 



those obtained by exchanging k with —k. 

A further simplification can be introduced by prop- 
erly choosing the distance between impurities L. For 



case the study of the state evolution greatly simplifies. 
Since at t — we had Iff) \ik)> at the time t we obtain 
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The correctness of the above approximation is checked 
by numerical integration of the Hamiltonian evolution. 
In Fig. \I\we compare the probabilities derived from the 
coefficients in Eq. <|12[1 with the exact results. The agree- 
ment between perturbation theory and numerical results 
is remarkable. 



IV. ENTANGLEMENT MEASURES 

The two-impurity reduced density matrix p is obtained 
by tracing out the electron degree of freedom: in the basis 
spanned by the states {|TT) , |iT) 7 iTi) j III)} we nave 
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Given p, we can compute the corresponding degree of 
entanglement by means of the concurrence G^ 2 . The 
concurrence between two qubits is defined to be C = 
max{Ai — A2 — A3 — A4, 0}, where A r is the square root 



of the rth eigenvalue of R 
matrix p is defined as p = 



-- pp in descending order. The 
(a y ® <Ty)p*((7y ® cry), where 
p* is the complex conjugate of p. Since the eigneval- 
ues of R are {0, 0, 0, (16/81) sin 4 (3 Jt/N)}, the corre- 
sponding A r are {0, 0, 0, (4/9) sin 2 (3 Jt/N)}. So we have 
C = 4/9 sin 2 (3 Jt/N). In Fig. |we report the numeri- 
cal value of C as a function of time and compare it with 
the analytical expression. In this way we have calculated 
the amount of entanglement arising spontaneously from 
scattering processes. Actually, there is a different strat- 
egy from which a higher degree of entanglement, the so- 
called localizable ent anglemen t 1 3 1 14 , could be extracted. 
It consists in a projective measurement performed on the 
electron degree of freedom. If we project ||T) \ik) t onto 
(life) + |T-fe)) /V2 we get the state 



(Tfe| + (T-fc| 
V2 
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As a result of the projective measurement, a maximally 
entangled state appears for all the times. However, the 
probability P of actually finding (\"\k) + lT-fc)) /\/2 is dif- 
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FIG. 2: Concurrence as a function of time. The solid line 
is derived from the theoretical model described in the text, 
while the dashed line corresponds to the numerical calcula- 
tion. The system parameters are the same as those defined in 
the caption of Fig. 
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That is, the process is probabilistic instead of deter- 
ministic. We obtain a success probability that evolves 
in time with the same law of C. The two results 
have the following interpretation. In the first case 
the two spins are spontaneously entangled by elec- 
tron scattering. The amount of entanglement is ob- 
viously related to the probability of finding the com- 
ponent CUT) |Tfe) + lit) |T-k) + III) ITfe) + |U) |T-fc»- In 



the physical procedure related to the projective measure- 
ment, we consider only this component. The same time 
evolution is then associated to two different kinds of pro- 
cesses, the first one being deterministic and the second 
one being probabilistic. In fact, without any kind of pro- 
jection, we have a limited degree of entanglement. By 
projecting on the electron state we can reach C = 1, but 
the price to pay consists in a limited success probability 
of the projection. The best strategy to adopt will depend 
on the specific application. 



V. CONCLUSIONS 

In this paper we have discussed the problem of entan- 
gling two distant spins embedded in a solid-state environ- 
ment through the interaction with a conduction electron. 
We have described explicitly what happens when a finite- 
size system is considered. By applying a weak-coupling 
approximation, that is, by neglecting nonresonant scat- 
tering states, we have solved analytically the evolution 
in time of the state associated to an incoming conduc- 
tion electron. In fact, the existence of discrete levels 
and the weakness of the coupling make possible a res- 
onance between the scattered states and one level of the 
energy band. A comparison between theoretical results 
and exact numerical results has been presented, showing 
the accuracy of the weak-coupling approximation. As a 
result of the evolution, two-spin entanglement appears. 
We have analyzed the emergence of entanglement merely 
from evolution through the reduced density matrix, and 
the role of projection on the electron spin state. In con- 
clusion, magnetic scattering in a discrete system has been 
shown to create entanglement both in a deterministic and 
in a probabilistic way. 
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